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This talk deals with two topics in multiparticle dynamics investigated by means 
of the DELPHI detector at LEP. Related to the first topic, we have used the 
data collected at 183 GeV to measure the average charged particle multiplicity in 
^ , e + e — — > bb events. The result is remarkably in agreement with QCD predictions, 

Q , while it is more than two standard deviations larger than calculations assuming 

i that the multiplicity accompanying the decay of a heavy quark is independent of 

1 the mass of the quark itself. The second topic deals with a recent theoretical model 

' by Andersson et al., in which soft gluons order themselves in the form of a helix 

I i at the end of the QCD cascades. In our data at the Z peak, we have found no 

evidence for such an effect. 

7— I ' 

> : 

t^J- , 1 Charged Particle Multiplicity in bb Events 

CO ; 

The study of the properties of the fragmentation of heavy quarks compared 
to light quarks offers insights in perturbative QCD. Particularly important is 
qq \ the difference in charged particle multiplicity between light quark pad. heavy 

■ miark initiated events in e + e~ annihilations. QCDjjredicts (see EoEru and 

""^ | El for a recent review) , somehow counter-intuitively □, that this difference is 

j^j ■ energy independent; this is motivated by mass effects on the hadronization. 

The present experimental teste, although preferring th^QCD-motivated 
scenario, were not conclusive (see tJ and references therein, Llfll§t2l). At LEP 2 
<-j ', energies, however, the difference between the QCD prediction and the model 

J" ■ ignoring mass effects is large, and the experimental measurement can firmly 

distinguish between the two hypotheses. 

Data corresponding to a total luminosity of 54 pb^ 1 collected by DELPHI 
at centre-of-mass energies around 183 GeV during 19117 were analysed. A 
description of-.the DELPHI detector can be found in £3; its performance is 
discussed in E3. 

Hadronic events were selected as in t3. The cross-section for e + e~ — > 
qq(7) above the Z peak is dominated by radiative qqy events, such that the 
effective hadronic centre-of-mass energy, y/p , is generally smaller lhan y/s. 
Only events with yfp (determined using the procedure described in Ii3) larger 
than 170 GeV were used. 

The influence of the detector on the analysis was studied with the full 
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DELPHI simulation program, —DELSIM III Events were generated with 
PYTHIA 5.7 and JETSET 7.4 M Partem Shower (PS). The particles were fol- 
lowed through the detailed geometry of DELPHI giving simulated digitisations 
in each subdetector. These data were processed with the same reconstruction 
and analysis programs as the real data. 

A 90% b enriched sample has then beep_obtained by requiring that the 
&-event tagging variable y defined as in Ref. lij was lower than 5 x 10~ 5 . From 
this sample the average bb multiplicity has been measured by unfolding via 
simulation for detector effects, selection criteria and initial state radiation. 
The value obtained was < n bb >= 29.54 ± 0.79(stat). _ 

The relative systematic uncertainty on < n bb > was assumed as in Q, so 
that: 

< n bb >= 29.54 ± 0.79(stat) ± 0.63{syst) . 

It was verified that the value above is stable within ± 0.13 units with respect 
to variations of the cut on the y variable. 

The difference 8u between the bb multiplicity and the multiplicity in 
generic light quark I — u,d,s events can then be extracted by using the 
equation: 

< rich >= Rb < n bb > +R C < n cc > +Ri < na > , 

where the fractions Rb = 0.16, R c = 0.26, Ri — 0.58 are the ones predicted by 
the Standard Model. For the cc multiplicity, the PYTHIA prediction of 27.6 
was assumed, with a systematic ermr of ±1.5 units. For < n c h >, we used 
the result published by DELPHI in H 

< n ch >= 26.58 ± 0.24(siai) ± OM(syst) . (1) 
Finally we obtain: 

6 b i = 4.23 ± 1.09(stat) ± lAO(syst) . (2) 

The value assumed for < n cS > is consistent within better than 1.5 standard 
deviations with both < > and < n bb >, neglecting the errors on < n i; - >, 

< n bb >: we thus verified a-posteriori that the systematic error assumed for 

< n cS > was conservative. The dependence of Sm on < n cS > is anyhow small, 
and the systematic error assumed gives a contribution of ±0.67 units to the 
systematic error on Sbi- 

Flavour-Independent Fragmentation — In a model in which the 
hadronization is independent of the mass, one can assume that the non-leading 
multiplicity in an event, i.e., the light quark multiplicity which accompanies 
the decay products of the primary hadrons, is governed by the effective energy 
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availabla-to the fragmentation system following the production of the primary 
hadrons One can thus write: 

S bl (E cm ) = 2 < n^ emv) >+J^ dx B f{x B ) J dx B f(x B ) n<r ((l- ^-±^ ) E cm ) 

~ n ll(Ecm) , 



where < n B ecay > is the average number of charged particles coming from the 
decay of a B hadron, xb (x b ) is the fraction of the beam energy taken by the B 
(B) hadron, and /(xb) is the b fragmentation function. 

We assumed 2 < n^ eca!/ ' >= 11.0 ±-,0.2 a, consistent with the average < 
measur pd at LEP 111 For f{x B ), we assumed a Peterson 
function of average 0.70 ±0.02 E3. The value of n t i(E) was computed from the fit to 
a QCD formula Ej including the resummation of leading (LLA) and next-to-leadiBg 
(NLLA) corrections, which reproduces well the measured chargedrSjiultiplicities E3, 
with appropriate corrections to remove the effect of heavy quarks lj. 

The prediction of the model in which the hadronization is independent of the 
mass is plotted in Figure ^ as a function of the centre-of-mass energy. There are sev- 
eral swiations of the model in the literature, leading to slightly different predictions 
(seeli3 and references therein). 

QCD Calculation — The large mass of the b quark, in comparison to the scale of 
the strong interaction, A ~ 0.2 GeV, results in a natural cut off for the emission of 
gluon bremsstrahlung, thereby avoiding the treatment of divergences in perturbation 
theory. Furthermore, where the centre-of-mass energy greatly exceeds the scale of 
the b quark mass, the inclusive spectrum of heavy quark production is expected 
to be well described by perturbative QCD in the Modified Leading Logarithmic 
Approximation (MLLA). | — . | — . 
The value of S b i has been calculated in perturbative QCDlilI: 

S b i = 2 < n ( g ecay) > ~ < n ir (e 1/2 m b ) > +0(a a (m b )) < n lt {m b )) > . 

Although the last term in the equation limits the accuracy in the calculation of 8 b i, 
one can conclude that S b i is approximately independent of E cm . 

The calculation of the actual value of 8 b i in u gives a value of 5.5 ± 0.8; aa 
approximation including the last term in the above equation gives a value of 3.68 u. 
This demonstrates the importance of the contribution 0(a s (m b )). A condition less 
restrictive is the calculation of upper limits. An upper limit 5u < 4.Viis given in u, 
while from phenomenological arguments, 8 b i < 4 is predicted in Ref. u. 

In Figure [j] the high energy prediction from QCD is taken from the average of 
the experimental values of 8 b i up to mz included, < S b i >= 2.96 ± 0.20. 

Our measurement of 5u is fully consistent with the hypothesis of energy in- 
dependence, and more than two standard deviations larger than predicted by the 
naive model presented in the beginning of this section. 
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Figure 1. The present measurement of 5;,; compared to previous measurements as a func- 
tion of the centre-of-mass energy, to the QCD prediction, and to the prediction from flavour- 
independent fragmentation. 



2 Is There Screwiness at the End of the QCD Cascades? 

Two colour charges moving apart from each other produce colour radiation in the 
form of gluons. At the end of the parton cascade one enters into a non-perturbative 
region, where perturbative QCD breaks up due to the small momentum transfer 
and hence large strong coupling constant a 3 . There are non-perturbative models 
to describe the transition regime between the perturbative region and final state 
hadrons, in which the largest contributing Feynman diagrams are chosen based on 
the coherence of the gluons. The gluon emissions are ordered by the transverse 
momentum in the dipole model EH, and by emission angle in the Webber-Marchesini 
model E3. The azimuthal angle of the gluons is taken to be uniformly distributed 
ignoring possible spin effects. 
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In a recent theoretical paper 13, it is studied what happens at the end of the 
QCD cascades: the soft gluons have a transverse momentum of the same order as 
the mass of the emitter, leading to a situation in which a large recoil may spoil the 
original coherence. With the assumptions that: 

• the effective coupling is large enough so that the number of emitted gluons is 
as large as possible; 

• the emissions conserve helicity; 

the Authors show that the closest packing of soft gluons is obtained by using the 
extra degree of freedom provided by the azimuthal angle and arranging the gluons 
along a helix. The helix is described by the parameter r = dy/dcf>, where y is 
rapidity and <f) is the azimuthal angle with respect to the original parton direction. 
The correlation between the azimuthal angle and rapidity cait-be tested by using a 
Fourier power spectrum, which is called screwiness S in Ref. Ej and defined as: 

S(lo) = £ P e | £ exp(i(a% - to))| a , (3) 

e j 

where the second sum goes over the gluons, the first sum goes over all the events, and 
P e is an event weight. Screwiness would then manifest itself as a peak in the power 
spectrum at w = 1/ <r>. If a large number of gluons were emitted isotropically, 
there should be no peak. Asymptotically screwiness should only depend of the 
number of emitted gluons at large values of uo. 

Experimentally one observes stable final state particles instead of gluons. More- 
over, in such an analysis it is convenient to use charged particles only, due to the 
better accuracy in the measurement of the momenta at low energy. Normalized 
screwiness can now be defined as 

s m=k — — -x;iE ex p«^-^-))i 2 , (4) 

J'eueiit \ J* particle 

e J 

where the second sum goes over the final state particles, and the first sum goes 
over all the events. Rapidity y and azimuthal angle 4> ar e defined with respect to 
the thrust axis. Ordering of the final state particles would produce a peak in the 
S(ui) spectrum, with a maximum peak size of S(u)max =< N par ti C ie > for maximal 
ordering. An isotropical particle emission would give asymptotically S(ui) ~ 1 at 
large values of u), and no peak. When u) approaches zero, screwiness measures the 
transverse momentum balance of the final state particles and it should approach 
zero. 

For a fixed \fs, a certain finite rapidity interval is available. The central rapid- 
ity plateau is not flat, producing fluctuations in the power spectrum. Selecting the 
central rapidity region and using only charged particles with experimental accep- 
tance cuts, the number of final state particles available to calculate the screwiness 
is reduced. 
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Figure 2. Power spectrum for simulated data samples with r = 0.7 (full line), r = 0.5 
(dashed line) and t = 0.3 (dotted line). Only stable charged particles with p > 0.15 GeV/c 
in \y\ < 2 were retained. 



A typical screwiness spectrum is shown in Figure |2j where the parameter r was 
set to 0.3, 0.5 and 0.7. Only charged particles with p > 0.15 GeV/c in \y\ < 2 were 
retained. A clear peak can be observed at uj = 1/ <r> when r = 0.7. With a 
r value of 0.5 there is still peak in the spectrum which could be observable, while 
with t = 0.3 the winding is too fast to be observed. The events were generated 
using a JETSET Monte Carlo simulation program, which did not contain parton 
showering, but in which a helix-ordering of the final state clusters was implemented. 
The average multiplicity of charged particles with the above selections was 5.2, which 
was roughly the same as with the default JETSET. 

With soft gluons included in the event generator, the transverse momentum 
fluctuations event-by-event cause that the screwiness does not decrease when u> goes 
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DELPHI 




Figure 3. Upper: power spectrum for DELPHI data (points) compared to JETSET with 
full detector simulation (histogram). The inner error bars on the data points correspond 
to the statistical errors, the outer ones to the sum in quadrature of the statistical and 
systematic errors. Lower: as for the upper part, for events with at least five charged 
particles and a missing momentum component transverse to the thrust axis smaller than 
1 GcV/c. 



to zero if one considers a narrow window in rapidity space with a few particles. The 
multiplicity of particles in the event affect the behaviour of the screwiness at small 
values of ui. At large values of ui the normalized screwiness approaches one. 

A sample of hadronic Z decays collected in 1994 were selected by requiring 
that the events contained at least five charged particles within the angular range of 
| cos6\ < 0.96, that the thrust axis was within the angular range of | cos#t| < 0.7, 
and that the thrust value was larger than 0.98. Charged particles were accepted if 
the measured momentum was greater than 0.15 GeV/c, track length greater than 
30 cm and the transverse and longitudinal impact parameters less than 5 and 10 
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cm, respectively. Particles within \y\ < 2 were retained for the calculation of the 
power spectrum. The final selected sample used for this analysis consisted of 50,000 
events. 

Screwiness of the selected events is shown in the upper part of Figure |^ together 
with the full simulation. The simulated data were generated with the JETSET 7.3 
parton shower simulation program, using input parameters reproducing the main 
features of the DELPHI data Ej. The generatsdi events were then simulated with 
the DELSIM full detector simulation program E3 and reconstructed as data. The 
data show no peak in the screwiness distribution, and the JETSET model, with 
random emission of soft gluons in the azimuthal angle, reproduces well the data 
distribution. 

The statistical error of the data was about 0.5%, computed as the standard 
deviation of a set of Monte Carlo simulation samples with the same population as 
the data. The systematic uncertainty was estimated to be around 1%. The main 
source of systematic uncertainty was expected to originate from the small differences 
in the rapidity distribution in data and simulation. The following procedure was 
used to estimate the inability of the simulated events to reproduce the data. The 
data were first weighted in such a way that the rapidity distribution was flat. The 
data were then weighted with the weights which would make the simulated rapidity 
distribution flat. The difference of the results, divided by the data distribution, was 
then the estimate of the relative systematic error. 

The same analysis was repeated by selecting only events with at least five 
charged particles with \y\ < 2, but no effect was seen. Finally, 10,000 events were 
selected with at least five charged particles with \y\ < 2 and such that the sum of the 
momenta of these particles had a momentum component transverse to the thrust 
axis smaller than 1 GeV/c (lower part of Figure ^|). The disagreement between data 
and simulation after these cuts is at the 5% level. The data does not show, however, 
a typical screwiness pattern. 

In conclusion, no evidence for screwiness was found, and results were consis- 
tent with random emission of gluons in the azimuthal angle. Given the DELPHI 
experimental sensitivity, it was found that the helix parameter r = dy/dcj) cannot 
be larger than about 0.5. 
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